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Mixed Path Ground Wave Propagation: 
1. Short Distances 

James R. Wait 

An expression is derived for the mutual impedance between two short vertical antennas 
on a flat earth with a straight boundary separating two media of differing electrical constants. 
After making some approximations that are valid at low^ and medium frequencies and where 
the antennas are not near the boundary, the integral formula for the field is evaluated for a 
wide range of the parameters. The numerical results computed in this paper are shown to be 
in reasonably good agreement with experiment. Finally, the effect of the obliqueness of the 
boundary is considered by a refinement of the stationary phase evaluation of the integrals. 

1. Introduction 

Coiisi(l(M-«il)h' inl(M"('st has Ixmmi sliowii rcccMilly in tlu^ propagation of ground waves over an 
inhomogoneous conducting earth. For niany purposes, such as estimating coverage of broad- 
cast transmitters, it is usually sufficient to assign an (^quivalcMit conductivity to the path. A 
simple and effective method for estimating this (HpiivahMit vahie from a conductivity profde of 
the path has recently been proposed by Suda [1].' Another technique described by Kirke [2], 
which is similar, has been called the (^qiiivahMit-distancc^ method and is also simple to a])pl3^. 
The best known method, however, is due to Millington |:^]. Although semiempirical it does 
appear to be valid for a wide range of freqiKMicies and ground constants. It also predicts the 
recovery efl'ect that occurs wIumi tlu^ wave ])asses over a boundary from an area of poor 
conductivity to one of good coiuhictivity. Millington (h)es not ()rovide a theoretical justification 
for his method, although he (k)es indicate^ that bis foiinulas for high frequencies are (compatible 
with the expected behavior of the height-gain functions over the media on both sides of the 
boundary. Clemmow [4] in an elegant (Ussertation obtains a I'igorous solution for a line source 
on flat earth parallel to a boundary sepai'ating two media. He makes a limitcul comparison of 
his formula with Millington^s method, and shows that the agreement is good. In a more 
recent theoretical approach to the subject, Bremmer [5] formulates the problem in terms of an 
integral ecpiation, which he solves by operational methods. After consideral)le nnuiipulation he 
succeeds in showing that his result is mathematically equivalent to that of CMennnow. 
Bremmer considers several limiting cases and also establishes the validity of Millington^s 
formidas at high frequencies. 

It is the purpose of the present paper to extend Bremmer 's result, with particular attention 
being paid to the phenomena at low and mediinn radiofrequencies where it does not seem 
possible to obtain convenient series expansions for the field by the Bi'emmer metliod. For 
the sake of completeness, the problem is reformidated in terms of mutual impedance between 
two antennas located on a flat earth, with a straight boundary separating two homogeneous 
media. In this case there is no doubt as to whether the reciprocity theorem is satisfied or not. 
Tlie line joining the two antennas makes some angle with the boimdary that should not be 
near zero. An integral equation is obtained foi' tlie field which is similar to one formulated 
by reinl)erg [6] for pi'0])agation over a rough groimd. In the present paper, the integral 
equation is solved by moncM'ical means for a range of parameters that are appropriate for short 
distances at low and mediiun frecjuencies. Jt is intended to extend the calculations in a later 
paper to situations where the earth's ciirv^atine must be considered; that is, for distances in 
miles greater than about 50X^'', where X is the wavelength in meters. 



' Figures in brackets indicate the literature rofdrences at the end of this paper. 
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2. Formulation 

The surface of the earth is considered to be flat and defined by z=0 in a simple cartesian 
coordinate system (x,y,z). Points abov^e the earth correspond to positiv^e values of 2. The 
earth mecUum to the left of a boundary line defined by y=x tan 6, as indicated in figure 1, has 
a conductiv^ity a and dielectric constant e. The earth medium for points to the right of the 
boundary line has a conductiv^ity of ai and dielectric constant ei. Short hertzian dipole antennas 
are located at the points A and B, which can be located anywhere along the x axis just abov^e 
the surface of the ground at 2=0. 

It is now assumed that a current_/a applied at the terminals of antenna A would produce 
electric and magnetic fields Ea and Ha for the case when the earth media were homogeneous 
with electrical constants a and e. A time factor ex'p{i^t) is implied. The mutual impedance 
between antennas A and B for this homogeneous case is denoted as Zab- When the earth becomes 
inhomogene_ous in the manner described above, the fields of antenna A with the same current 
la become Ea and //«, and the mutual impedance becomes Zat- It is convenient to regard the 
change from Zab to Zab as a consequence of the changes of the currents within the earth when 
it is transformed from a homogeneous to an inhomogeneous state. These currents will result 
in electric and magnetic fields Ea—Ea and Ha~Ha and a voltage laiZa^—Zab) at the terminals 
of antenna B. A current /& is now considered to be applied at terminals of antenna B, and 
the resulting electric and magnetic fields are E^ and H^ over the homogeneous earth of con- 
stants G and e. It then follows from Ballen tine's ^'corollary I" of the electromagnetic reci- 
procity theorem [7,8] that 

Za-Z,,= ,V f ( \EbXHa-E'aXHb]zdxdy, (1) 

J^aJ-h J J S 

where the integration extends over the whole ground plane S, and the subscript 2 indicates 
that the 2 or normal component of vector products is taken. 

A simplification is now made by introducing the concept of surface impedance [6,8], that 
is, the tangential electric and magnetic fields on the surface of the earth are assumed to be 
related by a complex constant of proportionality. More specifically, 

E^^^-vHA 

y (2) 

EyC^7]Hr J 

where rj, the surface impedance, is assumed to be equal to the value obtained for a plane wave 
at grazing incidence on the flat earth. It is given by 



ion '^^ 



!['+g"' 



where i8=27r/wavelength, and a=(i(T}j.co—eiJLco^y'-, with M = 47rX10~^ A subscript 1 is to be 
added to rj when the electric constants are ai and ei. 
Ecpiation (1) can now be written 

IJbAZ={rj- vOJj^^^^ {H'arHbt)dxdy, (4> 

where A.Z=Zab—Zab is the change of the mutual impedance from the situation of the homo- 
geneous earth of surface impedance -q to the inhomogeneous model indicated in figure 1. The 
quantity Hbt is the tangential magnetic field of the antenna B over the homogeneous earth, 
and Hat is the tangential magnetic field of the antenna A over the inhomogeneous earth. 




Figure 1. Schematic diagram of antennas A and B on a flat earth with a straight bound- 
ary separating two media with different conductivities. 



The integration now extends over the surface Sq to the rigiit of the boundary hne. Since 
H'at is not actually known, eq (4) is actually a surface integral equation formulation of the 
problem. Before attempting to solve it, a siini)lification can be made if tlve principle of sta- 
tionary phase is utilized. 

Letting the hertzian auteiuui at V> have an effective height K, the tangential magnetic 
field at P, a distance r from B, considering the earth as homogeneous and flat, is given by 






^^(l+i)i^Tr,r/)TixO, 



where F{r,rj) is Sommerf eld's ground-wave ath^nuation factor defined by 



(5) 



(6) 



where 



])-- 



'^ 2 a'\ ^a'J 2 ViSOtt/ 



is the numerical_ distance. Equation (6) is not exact, being an approximation valid for 
l3^<^\a\K Tr and '^^ are unit vectors in the directions of increasing r and z, respectively. The 
function F{r,7j) has been treated numerically by Norton [9] and others. The tangential mag- 
netic field Idat at a distance R from A is of the form 



,,;,= W« ,-.. (i +^) F'iB, ., .0 (-.Xl.), 



(7) 



where F'(R,ri,'r\^) is some function of 7?, r,. and t/j, and can be expected to be slowly varying 
compared to e~*^'*. is is a unit vector in the direction of increasing B. 



Using eq (5) and (7) and denoting the angle between the vectors v and Ir as 5, it follows that 
the expression for the mutual increment AZ can be written 

AZ=-(,-r;0 ^ J£ ^-^^ (l+^^) (l+^^) F{:r,,)XF^(r, ,,,,) cos My. (8) 



'5o 

^^here 



^==[(^0+^)^+2/^]^' and 

The major contribution to the integrand occurs when the phase of the exponential term is 
nearly constant, since the others are relatively slowly varying. With this in mind, the 
exponent is expanded in a powers series in if as follows: 

r^E^n+E^^y; [(-;-)+(^,)] for E.>.>0 



and 



where terms in y^, y^, etc., are neglected. A further approximation to the integrand in eq (8) 
is to replace (l + l/ip?-) and (l + l/i^R) by unity. This will be justified if the antenna A or B 
is not near the boundary (i. e., ^Eo and Pro':$>l). This is essentially equivalent to stating 
that the effect of the induction field of tlie antennas is not considered so far as the boundary 
is concerned. 

The integral expression for AZ now has the following approximate form: 



2i{ro+x){Ro-x)J 



The integration over x only extends from to Ro, because the contribution from the integrand 
for x>i?o is negligible due to the presence of the rapidly varying function 6~^'^^. 

3. Numerical Solution 

The mutual impedance Z between the two short antennas A and B for the homogeneous 
ground of electrical constants a and e is given by 



'27r(ro+i?o) 



It is now convenient to express the mutual impedance Z^ between the antennas on the inhomo- 
geneous ground as follows: 



where F' is the unknown attenuation function. The latter bracketed term can be replaced by 
unity because )S(ro+i?o) is large compared to 1. Employing eq (10), (11) and (12), it follows 
that 

F (, + i^,,,,,o^i^(/o+Ao,.)-ieo-(r;-,.;^-^-.^j J^ [(ro+x){Eo~x)Y^ ^^^^ 

Tliis is an integral equation for F\ It is immediately apparent when x<CO such that A and B 
are to the left of the boundary 

F'{ro+Bo,v,Vi)^F{ro+Bo,vy 

In otlier words, subject to our stationary phase approximation, the attenuation function F^ 
between A and B, when they are located on one side of the boundary, is characteristic only of 
the electrical constants a and e of the homogeneous ground between them. For exactly the 
same reason, the function J^'(ro+x,r;,r7i), where it occurs in the integrand of equation, for the 
general case of xy>0, can be replaced by F(ro-{-x,Tji), which is characteristic of propagation 
from A to the point ;/^O0) over a homogeneous ground of electrical constants ai and ei. 

After a cliarigo of vnriahl(\ tlic final (expression foi- F^ can be written in diniensionless form? 
yielding 



witli - 



and 



F{p)=l — i2py''e-^ I e-^'dot 
2 V1207r/ 



P«= 



K=^+}^^. and V= ^° 



-icx)e -/?o+'^o 

The above expression for F^ is then givcMi in terms of the attenuation function PXpo) , charac- 
teristic of propagation between A and B on a homogeiu^ous ground, and a correction term that 
accounts for the change of the electrical constants of the ground to the right of the boundary. 
It should be noted that po is the '^numerical distance^' between A and B with regard to the earth 
medium to the left of the boundary, whereas PoV/K is the numerical distance between the 
boundary to the point B with regard to the earth medium to the right of the boundary. Un- 
fortunately this expression, which involves an integration over products of error functions, is 
not readily expressible in closed form. However, the integral can be approximated in certain 
limiting cases. For example, when Pq{1 — V) and PoV/K are small compared to 1, the power 
series expansion for the functions F{p) and F(pQ—Kp) can be employed. This has the form 

F(p)=:l-i(7^py^-2p + i^^lTp^''± . (15) 

It then follows without difficulty that 



2 For large distances, the function F must be modified to account for earth curvature. This will be considered in section 2. 



plus terms in p^(^, pi, etc. Another special case is when 2>o is large compared to 1, so that only 
the leading term in the asymptotic expansion of F(pQ—Kp) need be considered; and fm-thermore, 
V is assumed to be small compared to unity. In this instance 

F(po~Kp) 1 1 



The power series formula for F{p) is then employed, yielding, after integrating term by term ; 
F' 



F^Po) 



2^(l-Vg) [u_i^|lfp,_2pTi^ 4 j)f_i^w p\_ 8 pi' 1 



(17) 



where pi=poVIK, 

There are probably other limiting cases which will enable the integration with respect 
to p to be effected; however, it is believed that for application to low and medium radiofre- 
quencies the values of po are in general neither large nor small compared to one. With this 
in mind, it is considered desirable to evaluate F^ by a numerical integration for a range of 
Pq, V, and K. It should be noted thai, for the general case, po and K are complex so there are 
actually five parameters to consider for the two-media problem. In this paper, attention 
will be restricted to frequencies where the displacement currents in the ground can be neglected. 
That is, the ratios eco/cr and eico/cri are assumed to be small compared to unity, therefore jh 
= 7r[(i?o4-/'o)A] (toco/o-) and K(=(Ti/a) are real quantities. This is usually justified for frequen- 
cies less than 1,000 kc for typical ground constants [101. 

Employing the numerical values of the function Fij)) and F{pq— Kp) . the integral in eq 
(14) is evaluated by a graphical method. The function F^ipQ, V,K) is then plotted as a func- 
tion of Po from 0.1 to 5 for various values of K and V in figures 2 to 11. It is believed that 
the results plotted in this form can be adapted to a large number of practical situations. The 
curves are not shown beyond Po^o, since this usually corresponds to higher frequencies where 
displacement currents are nonnegligible. 

It is of considerable interest at this stage to compare these numerical results with those 
computed using Millington's method [3]. His empirical formula for F^(po, V,K) in the nota- 
tion of the present paper would read 

F'(r, VK^ {F(r^lK\F(.n\lMkzYM^P$ZlETV n S^ 

F (l'o,F,if)=|_i^(i>„/i^)F(/>o)^^^^(j_^y^j^(p^^J • (18) 

Values computed from this formula for the cases K=2 and K=co are indicated by O ai^d 
#, respectively, on figures 2 to 11. The agreement is quite reasonable, and therefore further 
support is given to the validity of the Millington method, which has been known to predict 
correctly, within 1 db or so, the attenuation over certain mixed paths. It is believed, however, 
that differences are not of sufficiently small order to be neglected if fairly precise values of the 
fields in amplitude and phase are required. 



Cl 



u 

Q 

Z) 



CL 
< 



1.0 
.9 
.8 
.7 
.6 
.5 
A 
.3 
.2 
.1 
.0 



FIG. 2 


































/ 


































/ 


/ 


































/ 


/ 


































/ 


/ 




















V=0.0 










/ 


/ 


































A 


^ 


HASE 


: i 


) 






























/ 


Ph 








_^ 


















/\ 


/ 


/ 




















^ 


^ 








/ 


/ 






















AMPLITUDE |f' 


7^-. 




























^ 












\ 


\ 


















_^ 


/ 


y 


















N 














^^^i 




{C 






















N 














JR\ 


yE< 


5 IN 


DEF 


>EN 


DE 


N 


T 


OF K) 




N 


N 




















X 


X 








































^ 










































170° 
160° 
150° 
140° 
130° 
120° 
110° 
100° 
90° 
80° 
70° 
60° 
50° 
40° 
30° 
20° 
10° 
0° 
160° 
50° 
140° 
130° 
20° 
10° 
100° 
90° 
80° 
70° 
60° 
50° 
40° 
30° 
20° 
10° 
0° 



CL 

> 



en 

< 

X 

a. 



.3 .4 .5 .6 7 .8.9 1.0 
NUMERICAL DISTANCE p 



3. 4. 5. 



Figures 2 and 3. Attenuation function |F'|e'>' /or transtnission hetween A {or B) and B {or A) over a mixed 

path plotted as a function of the numerical distance po. 
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FiGUKES 4 AND 5. Attenuation function \F'\e''^' for transmission between A {or B) and B {or A) over a 

path plotted as a function of the numerical distance po. 
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Figures 6 and 7. Attenuation function \F'\eJf' for transtnission between A (or B) and B {or A) over a mixed 

path plotted as a function of the numerical distance Pq. 
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Figures 8 and 9. Attenuation function |/^'|e'>' for transmission between A (or B) and B {or A) over a mixed 

path plotted as a function of the numerical distance po- 
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Figures 10 and 11. Attenuation function {E'le'^'for transmission between A (or B) and B {or A) over a mixed 

path plotted as a function of the numerical distance Pq. 
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Figure 12. Field strength versus distance curves for various mixed paths. 

The conductivity in millimhos per meter to the left and right of the boundary Is shown on each curve. The experimental results were communicated to 

the author privately by P. A. Field, formerly of the Canadian Broadcasting Corporation. , From parametric curves; — - — -, computed by Millmg- 

ton method; O, experimental (P. A. Field); , curve, assuming homogeneous ground throughout {V=K'=0). 

4. An Application 

To illustrate the application of these parametric curves, the field strength in miUi volts per 
meter is shown plotted in figure 12 as a function of the distance (=ro-\~Ro) in miles for five 
different mixed paths. The curves are normalized so that the field strength at 1 mile is 
100 mv/m. Experimental values supplied by P. A. Field are also indicated. The curves 
calcidated by the Millington method are also shown in figure 12. The conductivities on the 
near and far sides of the boundary are indicated on the figures and are expressed in millimhos 
per meter. The agreement between the computed and experimental results is quite good. 
For these short distances the curves calculated by the Millington method fell slightly below 
the curves computed from the integral formula. In view of the uncertainty in the exact 
physical features of the ground, the difference between the two methods of calculation hardly 
seems significant for these situations. It is interesting, however, to note that the Millington 
method in four out of five cases underestimates the recovery effect to a greater degree than the 
integral method. The dotted curves in figure 12 correspond to the case where the ground is 
homogeneous throughout, and has the conductivity of the earth on the transmitter side of the 
boundary. The difference between the solid and dotted curves is a measure of the recover}^ 
effect, and it is apparent that the theory indicates a gradual transition at the boundary and 
does not show any sudden or transient features at the boundary. It would seem from the 
experimental results that there is actually some type of disturbance at the boundary between 
the two media which is not predicted by the theory. It is believed that a more rigorous 
evaluation of the integral equation is necessary to describe the nature of the field near the 
boundary, 

5. The Refraction Effect 

Although this problem, was formulated for the mutual impedance between antennas A and 
B on either side of an oblique boundary, the consequent path of stationary phase did not depend 
on the inclination angle 6. In other words, any refraction effects were neglected. It is the 
purpose of this section to revise the stationary phase evaluation of the integral to account 
for the changes of phase velocity between the two media. It is admitted that the procedure 
is not rigorous and is based on certain physical ideas that have their roots in geometrical optics. 

With regard to figure 1, the direct ray AOB between the antennas would represent the path 
of stationary phase if refraction effects were neglected. Of course, if the boundary was at a 
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right angle to AOB (i. e., ^=90°), it would be rigorously justified to assume AOB was the path 
of the stationary phase. It can be expected, however, that for the oblique boundary the 
path of stationary phase would be along a line APB, where P was displaced along the boundary 
from O by an amount q. Furthermore, it can be anticipated that q woidd be small compared 
to fo and i?o. The location of the point P can be best obtained by regarding q as a variable 
quantity in order to find when the total phase along APB is stationary or when it is a minimum. 
This is essentially a statement of Fermat's priiuuple in optics. The (hfi'erence between the 
phase along the patlis AOB and APB is now considered to be due to two factors: (1) the actual 
increase of path length, and (2) the change of the ratios of the path lengths for the two respec- 
tive media. For example, if tlu^ left-hand medium is more poorly conducting than the right- 
hand medium (i. e., i^>l), tlie phase velocity is relatively less to the left of the boundary, so 
that AP is less than AG. More explicitly, the total phase $ over the path APB can be written 

^(g)=/S(?';+/b+</>'(F,i^,/^o) (19) 

where 

r = [rg+g'-2ro^cos^]^ 

A=^[^?+<zH2A^or/cos^]^^ 



y= 



ro+i?o 



J A i^(r+Ko)/ V \ 



The function cI)(V,K,Pq) m the above expression is taken to he of the same form as the function 
<l>{V,K,po) previously computed. The general scheme is then to vary q to find wluM'e ^{(j) 
is stationary. For the present purpose, liowever, it is convenient to make some further 
approximations utilizing the fact that g is small. The phase fuiu'tion then becomes 



^{g)^fi(rARo) [l +^ (ff y^^^+<l>'iV,K, />„), (20) 



\vlier( 



and 



The terms containing higher powers in (q/d) have been neglected and, in fact, in the expression 
for po, the term containing (q/dy can be also dropped, since (l)^(V,K,po) is a relatively slowly 
varying function of po. It is now convenient to express </>'(F) as a Taylor expansion, as follows: 

Because q is small compared to ro+/?oj only the first term or two of the expansion are significant. 
The phase function now becomes 

$(g)^^(r„+i?o)+0'(F,i^j)„)-^^-^-|^^ 5(V)+^;^^^ x^ 0)', (22) 

where d(V) = —d<l>{V)ldV is essentially a positive quantity. Tlie right-hand side of the 
equation is a quadratic in q and has a minimum or stationary value when 

q_qo_ cose8iV)V(l~V)\ 

d d 2 sin' 9 w{ro+Bo) ' ^ 
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This value of g, namely, go, which makes $(g) stationary, is then a measure of the deviation of 
the direction of the phase of the ground wave as it crosses the boundary. The appropriate 
values of 5(F) can be obtained directly from the parametric curves of (t>^(po,VjK). Actually, 
it is a little more meaningful to consider the angles 8 and 8i, which are the angles OBP and OAP 
in figure 1. These can be called bearing errors, as they are a measure of the change of direction 
of the phase front relative to a direct unimpeded wave between A and B. In terms of goi the 
bearing errors, expressed in degrees, are given by 



and 



go 1 • .180 
- — r~w n — J^^ sin d — 



go 1 . .180 
di ^ — f^- Y^ sm 



Employing typical values of the parameters and letting K= co , values of 5 and di were com- 
puted and are: 
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It might seem surprising at first glance that these values are so small. This is as should 
be expected, however, because the absolute phase velocities of the ground wave over the two 
media differ only by a fraction of one percent [lOj. It is noted that the bearing errors are 
largest when the antenna is near the boundary, and where the separation between antennas 
A and B is 50 wavelengths or less. Of course, as the obliqueness of the boundary becomes 
more noticeable (i. e., ^<45°) the calculated values of 5 and 5i would be larger. One should 
be cautious, however, in attempting to apply this method to the case where the boundary 
makes a small angle with the line AB between the antennas. In such an instance the approxi- 
mation used in evaluating the basic integral by a stationary phase principle becomes invalid. 

6. Conclusion 

The results obtained in this paper provide a theoretical check of the Millington semi- 
empirical method of calculating the amplitude and phase of a wave crossing a boundary sepa- 
rating two media. The parametric curves presented here should be convenient for making 
predictions for propagation along a two-media path at low and medium frequencies. It is 
also indicated that the refraction effects for a wave crossing the boundary obliquely are very 
small. 

I record my thanks to Loris Perry, who carried out most of the computations. 
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